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Abstract 



The Bose-Einstein condensates of alkali atomic gases are spinor fields with 
local "spin-gauge" symmetry. This symmetry is manifested by a superfluid 
velocity u s (or gauge field) generated by the Berry phase of the spin field. In 
"static" traps, u s splits the degeneracy of the harmonic energy levels, breaks 
the inversion symmetry of the vortex nucleation frequency fi c i , and can lead 
to vortex ground states. The inversion symmetry of O c i, however, is not 
broken in "dynamic" traps. Rotations of the atom cloud can be generated by 
adiabatic effects without physically rotating the entire trap. 



The recent discoveries of Bose-Einstein condensation in atomic gases of 87 Rb 7 Li ||, 
and 23 Na have achieved a long sought goal in atomic physics. They have also provided 
condensed matter physicists opportunities to study interacting Bose systems at a wide range 
of densities. The realizations of these condensates are made possible by the invention of a 
number of special magnetic traps, which trap atoms with hyperfme spin (F = 2) maximally 
aligned with the local magnetic field B. The reported Bose-Einstein condensations (T| |2j 
H are found in these (adiabatic) spin states. 

An immediate question is whether these alkali condensates differ from the familar 4 He 
condensate in any fundamental way. Unlike the spinless 4 He atoms, the trapped alkali atoms 
are in the F = 2 hyperfme spin state. Their condensates are therefore spinors of the form 




(1) 



where ipm is the field operator, m is a label for F z , (—2 < m < 2), $ is a scalar, and is a 
normalized spinor. Since the hyperfine spins are aligned with the magnetic field, ( is given by 
B • F( = 2(, where F is the hyperfine spin operator. The dynamics of < ip m > is therefore 
completely specified by that of the scalar field $, as in 4 He. One might then conclude 
that apart from extrinsic factors like density and external potential, there is no intrinsic 
symmetry difference between 4 He and alkali condensates. This is in fact the starting point 
of all current theories, which model the alkali systems as interacting dilute spinless Bose 
gases in harmonic potentials. Within these models, the effective Hamiltonian for the scalar 
$ has a global U(l) gauge symmetry, as in 4 Be. 

The actual symmetry of the spinor field (eq.(|l|)), however, is much larger than U(l). 
We call it local spin-gauge symmetry. It represents that a gauge change e lx<yyL,t ^ of < ip m > 
can be undone by a local spin rotation e~^ X//F - )B ( x '*' ) ' F . Because of this symmetry, the exact 
form of the effective Hamiltonian of the scalar <3> is not that of 4 He, but that of a neutral 
superfluid in a velocity field u s , or an electron in a vector potential potential A. The velocity 
(or gauge field) u s arises from the Berry phase of the spin field (. It is a direct reflection of 
the underlying spin-gauge symmetry. The purpose of this paper is to discuss various forms 
of spin-gauge effects. 

As we shall see, the strength of the gauge field u s is proportional to the gradient of 
the magnetic field B. It is generally quite small because B is generally fairly uniform in 
the trapping region. While these may be a justification of treating the alkali system like 
dilute 4 He, we note that spin-gauge effects can be magnified rapidly by variations of trap 
parameters and particle numbers. As we shall see, despite its weakness, the effect of u s can 
be observed over a wide range of trap parameters in both normal and superfluid phases. 
Moreover, in the condensed phase, particle interactions can reduce the excitation energy 
so much that the spin-gauge effect are significantly magnified. In particular, when the 
excitation energy lies below that of spin-gauge effect, the ground state will change abruptly 
to one reflecting this symmetry. Generally, spin-gauge effect increases as one moves away 
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from the center of the atom cloud. Understanding of this effect is therefore important in 
view of the current effort to produce larger and larger clouds. It is also important to note 
that nonuniformity of the magnetic field, however small, is what leads to magnetic trapping. 
Spin-gauge effects are therefore intrinsic properties of magnetically trapped atomic gases. 
Our major findings are : 

I. The adiabatic spin field ( of the alkali system generates a superfluid velocity u s . Its 
vorticity fl s = |V x u s is specified by the magnetic field through a topological term, (eq.(||) 
below). 

II. For cylindrical traps with static fields (or simply "static traps"), u s generates a Coriolis 
force which splits the degenerate harmonic energy levels. This splitting is independent of 
particle interaction and can be observed in both normal and superfluid phases. 

III. For static traps, the "background" rotation — fl s breaks the inversion symmetry of the 
vortex nucleation frequency O c i. For appropriate trap parameters, the system has a vortex 
ground state even in absence of external rotation. 

IV. For cylindrical traps with time dependent fields, (or simply "dynamic traps"), Q s has 
a quadrupolar structure. As a result, the inversion symmetry of O cl is restored. 

V. Adiabatic effects furnish a simple means to generate rotation of the atom cloud without 
physically rotating the entire trap. 

To begin, we first discuss the effective Hamiltonian. For brevity, we shall call hyperfine 
spins "spins" . The Hamiltonians of the alkali systems are of the form H = H s + V, where 
H s = J cfx?/>+(x) [— ^V 2 — /i a B(x, t) ■ F] mn ^ n (x) is the single particle Hamiltonian, M and 
A*a = — A*b/2 are the mass and magnetic moment of the atom, \ir is the Bohr magneton, and 
the factor 1/2 is the g-factor of the alkali atom. B is a sum of magnetic field configurations 
which can be static or dynamic V is the two-particle interaction between the atoms. To 
form a trap, the Zeeman energy —fi a B (or their time average) must behave like a potential 
well. If {C^} are the spin eigenstates along B, (B • = nC )> (—2 < n < 2), then the 
Zeeman energy — /i a B • F reduces Z7^(x, t) = |n/xsS(x, t) for the states C^. If hbB is an 
attractive well, U^ n ' is confining (deconfining) for n > (n < 0). This means that spin-flips 
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between n > and n < states can cause atoms to leave the trap. Since V generally cause 
spin flips unless both atoms are in the maximum spin state along the same quantization axis, 
(in which case spin flips are prohibited by angular momentum conservation), it depletes all 
but the "adiabatic" spin states in the trap. The resulting system is an interacting Bose 
gas with spins aligned with the local field B(x, t). 

To construct a theory for the adiabatic spin states, we expand ip m in terms of the 
spin eigenstates (( n \ tp m ('x,t)= Y,n=-2(m\ x i <P n ^( x >0- Expressing B = zcos/3 + 
sin/3 (xcosa + ysina), the explicit form of is 

=< m \U\n >, U = e -i«F Ze -if3F ye -i X F z ( 2 ) 

where F z \n >= n\n >. % is arbitrary. It is the gauge degree of freedom of the system, and 
is usually chosen to make the spinor single valued. The effective Hamiltonian TC can 
be obtained by rewriting the equation of motion ihd t x(j m = [tjj m , H] in the form ihd t (p^ = 
[4>( n \Tl]. One then finds H = H a d + H na d + Ti-etc- 'Had, referred to as the "adiabatic" 
Hamiltonian, contains cj)^ only. Tl naa - is the spin-flip (or nonadiabatic) Hamiltonian which 
consists of cross terms between (j)^ and <fi( n ^ 2 \ H e t c describes the transitions between 
different n ^ 2 states and can be ignored. Denoting cj)^ and as (j) and C, respectively, 
we have 

1 /fiV \ 2 

where U= U^= fi B B(x), W= (h 2 /2M)[\V(\ 2 + (C + V() 2 ] -ih(+d t (, V is the projection of 
V onto the adiabatic spin states. It is of the form V— JV(x — y)0 + (x) + (y)0(y)0(x), 
where V(x — y) is a short range potential. The velocity u s is defined as 

Mu s = -C+VC (4) 

Eq.(|3]) describes a Bose fluid in a "background" velocity field — u s , or a charge e system in 



Had = I dx< 



V, (3) 



a vector potential A if Mu s = eA/c. Under a local spin rotation exp iQ ■ F%(x 



Us 



u s + (F/i/M)Vx(x), which is equivalent to a local gauge transformation 0— » exp(iFx(x))0. 



This is a reflection of the underlying spin-gauge symmetry of (. The integral J c u s ■ ds is 
the Berry's phase of ( around a loop C. It can be easily calculated from the vorticity (O s ) 
of u s , which satisfies the Mermin-Ho relation [|J, 

n s = ^Vxu s = t aM B a VB p x V5 7 . (5) 

Eq.([5|) shows that the spatial variations of B necessary to produce the trapping potential 
will inevitably generate to a non- vanishing superfluid velocity u s . 

In the rest of this paper, we shall focus on the phenomena associated with the adiabatic 
spin fields, (described by TC a d only). Nonadiabatic effects will be discussed elsewhere 0. 
Typically, nonadiabatic effects of the trap can be ignored if its "Dirac center" is sufficiently 
far away from the atom cloud. The "Dirac center" is the point where B = and that 
the unit vectors B surrounding D wraps around the unit sphere n times, (n is a nonzero 
integer). If D resides in the cloud, the adiabatic spin field around D will develop a line 
singularity emerging from D, (a Dirac string), which will cause a lot of spin-flips. As we 
shall see, increasing field gradients enhances spin-gauge effects but at the same time moves 
D closer to the cloud. The field parameters discussed below are all within the range to keep 
the Dirac center sufficiently far away from the cloud. 

To be concrete, we consider "static traps" of the form, (V • B = V x B = 0), 

B(x) = B£ + G 1 (xk-yy) + Q- 

where r = (x,y), G\ and Gi are the first and second order field gradients respectively. 
Magnetic trap of the form eq.(||) is similar to that used in the 7 Li experiment 0. It is 
convenient to express the field gradients as G\ = B (j/L), G2 = B /L 2 . The trapping 
potential IA in eq.(|J) can then be expressed as 

+ 0|x/L| 4 . (7) 

where TiVL Zee = fx B B , u 2 z = fi B B /(ML 2 ), u z /lu ± = {j 2 - 1/2)" 1 / 2 = A, 7 2 > 1/2. For later 
use, we denote the longitudinal and transverse width of the ground state Gaussian of the 
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u = m Zee + -M 



2 2 1 22 



harmonic well (eq.(0)) as a z and a±, where a z = (h/Muz) 1 / 2 , a± = (h/Muj^) 1 / 2 . Typically, 
a z ,a ± « L. It is straightforward to show that W= (Ti 2 /2M)([sin/3V 'a] 2 + [V/3] 2 ), where 
a, (3 are polar angles of B as defined earlier. This term is smaller than the harmonic potential 
in U by a factor ( / -fa±/L) 4 and can be ignored in general. 

^From eq.(||), it is straightforward to show that [with \ — ~ a i n e Q-dl)]? u s= 
f (1 - B z /B) V[t a n-\B y /B x )], and 

u s = -^g) 2 zxr + 0(x 2 /L 3 ), n s = -zA(l) 2 + (|x|/L 3 ). (8) 

Thus, for |x| < L, spin-gauge effect generates a constant effective "rotation" — f2 s along z. 

An immediate consequence of Q s is that it generates a Coriolis force on the alkali system. 
This force can be detected by the applying an a.c. magnetic field along x, b = 6e~ iw *x. This 
field will generate a term (/j,Bb/2L)xe~ luJt in the effective Hamiltonian 7i a di as if a time 
dependent force f = (/i j e&/2L)e _ * a; 'x is present. It is easy to see that the equation of motion 
of the center of mass in the xy-plane, R = / + r0, r = (x,y), assumes the form 

M—f^- = —MuVR + 2M^ x a + f • (9) 

dt 2 ± dt K J 

which has resonances at to = uj± ± Q s (for uj± » Q s ). The degenerate clockwise and 
counterclockwise harmonic modes uj± are split by the Coriolis force. This splitting exists in 
both normal and superfluid phases, and can be easily shown to be independent of particle 
interactions. 

More pronounced effects can be found in the superfluid phase of alkali atoms with positive 
scattering length a > 0. Because of spin-gauge symmetry, the ground state energy functional 
becomes 



£($) ' 



2M 



— + Mu s J $ 



2 + (W + W)|$| 2 + ^|$| 4 . (10) 



When u s is small, eq.(0) can be written as £ u s )= £($, 0)-VL s C z , (C z = -i$*z- r x V$), 
which is the Hamiltonian density of a scalar superfluid in a container rotating with frequency 
Q s z. Let f2°L denote the vortex nucleation frequency in the absence of spin-gauge effect (i.e. 
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Q s = 0). Because of the "background" rotation Q s , the actual vortex nucleation frequencies 
ttfi for vortices with 2n circulation around ±z will be = VL° cl =p VL S . In particular, when 
Q s > hence < 0, vortex ground state will emerge in the absence of external rotation. 

The value of VL° cl has been studied for harmonic traps by a number of authors || [[7]. 
Using Thomas- Fermi approximation (TFA), which is good at large N ||, Baym and Pethick 
have shown that Vt° cl is reduced by particle interactions from its non-interacting value u± as 

n° cl /uj ± = Q- 2 \nQ 2 , Q = R±/a± = (15XNa/a±) 1/5 . (11) 

where R± is the transverse width of the condensate. Since Q° cl oc N~ 2 '\ Q s oc N°, the 
inversion asymmetry of the nucleation frequencies, (fl^ 1 —fl~ l )/(Q^ l +Q~ l ) ^flg/fl^, increases 
as N 2 / 5 . Thus, for sufficiently large N, the condition of vortex ground state f^/O^ > 1 can 
always be met. From eq.([8D and eq.(^), one finds that the ratio Qg/fl^ increases as the 
externally controllable parameters N, G\,B~ l increase. 

Figure 1 shows the ratio VL s /VL° cl as calculated from eq.([8]) and ( [TT]) for 23 Na, which has 
a positive scattering length a = 4.9nm. The asymmetry of the trap is set at A = 1/2. 
Four cases are considered : N = 5 x 10 5 (broken line) and N = 5 x 10 6 (solid line); and 
B Q = 3 and 5 Gauss, (denoted the numbers 3 and 5 respectively). They are chosen to 
indicate the direction of increasing spin-gauge effect as well as the conditions for vortex 
ground states. The field gradients considered extend to the Tesla/cm range. Although 100 
times higher than those in current experiments (~ lOOGauss/cm), they are easily achievable 
using superconducting magnets. At present, the largest iV produced is 10 5 . However, since 
iV has increased from 10 3 to 10 5 in last seven months, it is conceivable that condensates 
with iV = 10 6 can be realized in the near future. 

^From figure 1, we see that for N = 5 x 10 6 , asymmetry (flg/fl^) up to 10% already 
occurs around 6500 Gauss/cm, and reaches 1 (i.e. vortex ground state) around 5 Tesla/cm. 
It should be noted, however, eq.(||), (|7|), ([□]) neglected terms higher order in |x|/L, which 
is justified only when R z ~ 2R± < L, where R z is the longitudinal width of the condensate. 
For sufficiently large G\ or N, this condition will fail, at which point higher order terms in U 



and u s , as well as W begin to contribute and that the system will lose cylindrical symmetry. 
The values of G\ at which R z = L are marked by circles on the curves in figure 1, indicating 



that eq. (|TT|) is only accurate to the left of the circle. 

To calculate fi cl in the regime R z > L, we have calculated the energies (E± and Eq) of 
a ±27r vortex and the no vortex ground state using the full expressions of U, W, and u s . 
The nucleation frequencies Qf x are related to these energies as ~ (E± — E )/h ||. Our 
calculation are performed within TEA, which in the present context amounts to replacing the 
kinetic energy by {h 2 /2M){V9 + Mu s /M) 2 |$| 2 and ignoring the (ft 2 /2M)V|$| 2 term. We 
have minimized this approximated energy subjected to the constraint of constant particle 
number N = J |$| 2 , where |$| is the magitude of the order parameter of a ±2tt vortex (and 
the no vortex state) in the E± (and E ) calculation. Our results are shown in Figs. 2a and 
2b, which plot the ratio rf 1 = Q^i/uj± as a function of G\. We see that vortex ground states 
(i.e. r] + = 0) emerge around 5 Tesla for N = 5x 10 6 . Even though the other cases considered 
have not reached vortex ground state, they show strong broken inversion symmetry in the 
nucleation frequency, i.e. (r] + — f]^)/(f] + + T]~) are close to or over 50%. Finally, we note 
that the ratio Q s /ujj_, which describes the amount of energy level splitting in eq.([]) is of the 
order of 10~ 3 for the range of parameters we considered, (which is easily verified from the 
expression of a± and Q s ). This splitting, though small, is within the limit of detectability. 
Note also that the Dirac centers of the static trap eq.(^) are located at (±La/87 2 + 2, 0, 27L), 
(0, ±Lv / 87 2 + 2, — 2jL). For the cases we considered, A = 1/2, hence 7 = V4.5, these Dirac 
centers are quite far away from the center of the cloud as R z ~ 2R± is less than 2L for all 
cases considered. 

To further illustrate the spin-gauge effect, we consider dynamic traps like those in the 
Rb experiment 

B(x,t) =B [n(t)+b(x)], b(x) = (r - 2zz) /L, (12) 

where B /L is the field gradient, and n(t) = pcosu; t + qsinci; t is a unit vector rotating in a 
plane perpendicular to 1, and (p, q, 1) form an orthogonal triad. The adiabatic Hamiltonian 
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7i a d is now periodic in time with frequency u . Expanding 7i a d{t) in Fourier series of e mu) ° t i 
the time averaged (i.e. n — 0) term give rise to a static Hamiltonian of the form eq.(0) with 
U replaced by U = (i B B(l + [b 2 + (b ■ 1) 2 ]/(4L 2 ) + ..). When 1 = z, U has the cylindrical 
symmetric form eq.(^) with u\ = ^bB /2ML 2 , (uj z /uj±) 2 = 8. The rotational frequency 
uj has to fall in the range uj z < uj < ttzee for the spins to follow the local magnetic field 
adiabatically [|]. Using eq.@, the time averaged of the velocity field u s and Q, s associated 
with eq.(^) are found to be 



u s = 



yML*;™*'- n - = ims) {r - 2 ^- (13) 

Unlike the Q s of the static trap, eq.(|J), Q s is a quadrupolar field which has mirror symmetry 
about the xy-plane. As a result, = f^, which is similar to 4 He but is entirely different 



from the static trap |10| . To verify this symmetry in the regime where f2 s has similar strength 
as fl s in the static trap, the system has to be rotated up to the critical frequency Q c \, which 
can be as high as 100 rad/sec for G = 5000 Gauss/cm, B = 5 Gauss, and iV = 10 5 . While 
it is impractical to rotate the entire trap at such high frequencies, we point out below a 
simple way to rotate the trapping potential using the adiabatic effects. 

Consider the case where 1 deviates slightly from z, hence causing U to deviate slightly 
from cylindrical symmetry. If 1 precesses about z with frequency u p , u p « uj_ < u , U 
will rotate about z with the same frequency. (The time average now is to be understood 
as averaging over times faster than 2ir/u p ). Since hi deviates only slightly from cylindrical 
symmetry, the corresponding vector fields u s and Q s are essentially given by eq . (|13D . Since 
the rotation (u l) and the precession (u p z) of n can be generated by electromagnetic means 
from sets of stationary coils, rotation of the trapping potential can therefore be generated 
without physically rotating the entire trap. In a similar fashion, the potential U of the static 
trap can be made to rotate about z by the application of a small magnetic field rotating in 



the xy plane with frequency uo p « uo± ||11|| . Rotating the trapping potential in this fashion 
allows one to study the inversion asymmetry of i n both low and high field gradient 
regime. 



We have thus established Statement I to V. Our discussions also show that spin-gauge 
effects assume different forms in different traps. It is therefore conceivable that they can 
be made more prominent at lower field gradients by other ingenious design of traps. In a 
broader sense, the spin-gauge effect is only a subset of a much larger class of phenomena 
associated with the topological excitations of the spin field, which are suppressed by the 
magnetic field in the current traps. Should it be possible to release part (or all) of the spins 
degrees of freedom in a new trapping design, the spin-gauge phenomena of the resulting 
condensate will be truly remarkable indeed. 

TLH would like to thank Greg Lafyatis for discussions. This work is supported in part 
by NSF Grant No. DMR-9406936. 
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Caption 



Figure 1. The ratio £l 8 /£l cl as a function of field gradient G 1 for 23 Na in a trap with 
asymmetry A = 1/2. The solid and broken lines represent particle numbers N = 5 x 10 6 
and N = 5 x 10 5 respectively. The labels 3 and 5 on the figure denote B Q = 3 and 5 Gauss 
respectively. On each curve, the region to the left (right) of the circle indicates the condition 
R z < L{> L). 

Figure 2a (2b) shows the ratio rf 1 = Qf 1 /uj± a function of field gradient G\ for B Q = 3 
(5) Gauss. The meaning of solid and broken lines as well as the symbols "3" and "5" are 
identical to those in figure 1. 
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